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Abstract 



o 

(-H ■ The problem of finding "small" sets that meet every straight-line which intersects a given 

C^ I convex region was initiated by Mazurkiewicz in 1916. We call such a set an opaque set or 

a barrier for that region. We consider the problem of computing the shortest barrier for a 

CO \ given convex polygon with n vertices. No exact algorithm is currently known even for the 

simplest instances such as a square or an equilateral triangle. For general barriers, we present 

an approximation algorithm with ratio i + ^-'^^ = 1.5867 . . .. For connected barriers we achieve 

Cn ■ the approximation ratio 1.5716, while for single-arc barriers we achieve the approximation ratio 

r\ , ^i| ~ 1.5834. . .. All three algorithms run in 0{n) time. We also show that if the barrier is 

restricted to the (interior and the boundary of the) input polygon, then the problem admits a 

^ , fully polynomial-time approximation scheme for the connected case and a quadratic-time exact 

algorithm for the single-arc case. 

^ ■ Keywords: Opaque set, opaque polygon problem, point goalie problem, traveling salesman 

00 I problem, approximation algorithm, Cauchy's surface area formula. 

(N. 

^. ! 1 Introduction 

in. 

^5 ■ The problem of finding small sets that block every line passing through a unit square vi^as first 

considered by Mazurkiewicz in 1916 [3l]; see also O |22]. Let C be a convex body in the plane. 
Following Bagemihl j3], we call a set B an opaque set or a harrier for C, if it meets all lines that 
intersect C. A barrier may consist of one or more rectifiable arcs. It does not need to be connected 

r> ! and its portions may lie anywhere in the plane, including the exterior of C; see [3j, [7]. 

d \ What is the length of the shortest harrier for a given convex body C ? In spite of considerable ef- 

forts, the answer to this question is not known even for the simplest instances of C, such as a square, 
a disk, or an equilateral triangle; see t8j, [9, Problem A30], [13], [15], [16], [191 Section 8.11], [23] 
Problem 12]. The three-dimensional analogue of this problein was raised by Martin Gardner [20] : 
see also [21 [7]. Some entertaining variants of the problem appeared in different forms \2.b\ |29| 130] . 
for instance: What should a swimmer at sea do in a thick fog if he knows that he is within a mile of 
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a straight shoreline? The shortest known solution resembles the shortest known single-arc barrier 
for a disk of radius one mile; see [9l Problem A30]. 

A barrier blocks any line of sight across the region C or detects any ray that passes through 
it. Motivated by potential applications in guarding and surveillance, the problem of short barriers 
has been studied by several research communities. Recently, it circulated in internal publications 
at the Lawrence Livermore National Laboratory |10j . The shortest barrier known for the square, 
of length 2.639 . . ., is illustrated in Fig. [Upright). It is conjectured to be optimal. The current best 
lower bound is 2, established by Jones p4] . 
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Figure 1: Four barriers for the unit square. From left to right: 1: single-arc; 2-3: connected; 4: disconnected. 
The first three from the left have lengths 3, 2^2 = 2.8284 . . ., and 1 + \/3 = 2.7320 . . .. Right: The diagonal 
segment [(1/2, 1/2), (1, 1)] together with three segments connecting the corners (0,1), (0,0), (1,0) to the 
point (5 - #, ^ - #) yield a barrier of length ^/2 + ^ = 2.639 . . .. 



Related work. The type of curve barriers considered may vary: the most restricted are barriers 
made from single continuous arcs, then connected barriers, and lastly, arbitrary (possibly discon- 
nected) barriers. For the unit square, the shortest known in these three categories have lengths 3, 
1 + \/3 = 2.7320 . . . and ^/2 + ^ = 2.6389 . . ., respectively. They are depicted in Fig. E Inter- 
estingly, it has been shown by Kawohl [26] that the barrier in Fig. (TJ^right) is optimal in the class 
of curves with at most two components (there seems to be an additional implicit assumption that 
the barrier is restricted to the interior of the square). For the unit disk, the shortest known barrier 
consists of three arcs. See also [HI [19] . 

If instead of curve barriers, we want to find discrete barriers consisting of as few points as 
possible with the property that every line intersecting C gets closer than e > to at least one of 
them in some fixed norm, we arrive at a problem raised by Laszlo Fejes Toth |171I18| . The problem 
has been later coined suggestively as the "point goalie problem" [40] . For instance, if C is an axis- 
parallel unit square, and we consider the maximum norm,, the problem was studied by Barany and 
Fiiredi [3], Kern and Wanka ^8j, Valtr |43j, and Richardson and Shepp [40J. Makai and Pach |33j 
considered another variant of the question, in which we have a larger class of functions to block. 

The problem of short barriers has attracted many other researchers and has been studied at 
length; see also [U [T^ [23l [3T]l32j. Obtaining lower bounds for many of these problems appears to 
be notoriously hard. For instance in the point goalie problem for the unit disk (with the Euclidean 
norm), while the trivial lower bound is 1/e, as given by the opaqueness condition in any one 
direction, the best lower bound known is only 1.001/e as established in |40j via a complicated 
proof. 



Our results. Even though we have so little control on the shape or length of optimal barriers, 
for any convex polygon, barriers whose lengths are somewhat longer can computed efficiently. Let 
P be a given convex polygon with n vertices. 



1. A (possibly disconnected) segment barrier for P, whose length is at most ^+ ^^ = 1.5867 . . . 
times the optimal, can be computed in 0{n) time. 

2. A connected polygonal barrier whose length is at most 1.5716 times the optimal can be 
computed in 0{n) time. 

3. A single-arc polygonal barrier whose length is at most ^i| = 1.5834. . . times the optimal 
can be computed in 0{n) time. 

4. For interior single-arc barriers we present an algorithm that finds an optimal barrier in 0{'n?) 
time. 

5. For interior connected barriers we present an algorithm that finds a barrier whose length is 
at most (1 + e) times the optimal in polynomial time. 

It might be worth mentioning to avoid any confusion: the approximation ratios are for each 
barrier class, that is, the length of the barrier computed is compared to the optimal length in the 
corresponding class; and of course these optimal lengths might differ. For instance the connected 
barrier computed by the approximation algorithm with ratio 1.5716 is not necessarily shorter than 
the (possibly disconnected) barrier computed by the approximation algorithm with the larger ratio 

1 _^ 2+v^ ^ 1.5867.... 

However, we believe that the approximation ratios of the first two algorithms mentioned above 
are substantially better than 1.57. In support of this belief, we present a couple of lower bound 
examples for which the ratios are below 1.1. 

2 Preliminaries 

Definitions and notations. For a polygonal curve 7, let I7I denote the length (or weight) of 7. 
Similarly, if F is a set of polygonal curves, let |F| denote the total length of the curves in F. As 
usual, when there is no danger of confusion, we also denote by |^| the cardinality of a set A. We 
call a barrier consisting of segments (or polygonal lines) a segment harrier. In order to be able to 
speak of the length i{B) of a barrier B, we restrict our attention to barriers that can be obtained 
as the union of finitely many rectifiable curves. We first show (Lemma [1]) that the shortest segment 
barrier is not much longer than the shortest rectifiable one. 

Lemma 1. Let B be a barrier of length i{B) < 00 for a convex body C in the plane. Then, for any 
e > 0, there exists a segment barrier B,, for C, consisting of finitely many straight-line segments, 
such that £{Be) < i{B) + e. 

Proof. Since every line blocked by a rectifiable curve must intersect the convex hull of the curve, and 
since any rectifiable curve 7 : [0, 1] — )• M^ is a barrier for its convex hull conv(7) = conv(7([0, 1])), 
it suffices to prove the lemma for barriers consisting of a single simple (not self-intersecting) curve 
7, with C = conv(7). 

We can assume without loss of generality that the endpoints 7(0) and 7(1) of 7 belong to 
Bd(C), the boundary of C. Otherwise, we can shorten 7 by clipping its initial and final pieces 
before its first and after its last point that belongs to Bd(C). 

A minimal interval (by containment) [to, ti] Q [0, 1] with the property that 7(^0), 7(^1) £ Bd(C) 
and the points 7(to) and 7(^1) separate 7(0) from 7(1) along the boundary of C is called a shifting 
interval. Obviously, no two shifting intervals share an interior point. For each shifting interval 
[to, ti], replace the arc 7([to, ti]) by the straight-line segment 7(to)7(ti), to obtain a rectifiable curve 



7o : [0, 1] — )• K^ with ^(70) < ^(7)) which is a barrier of C. The segment 7o([iO)ii]) = 7(^0)7(^1) is 
said to be a chord. We also have conv(7o) = conv(7) = C. 

The curve 70 is composed of arcs that are alternately chords and (possibly one-point) pieces of 
the boundary of C. That is, there is a (possibly two-way infinite) sequence of numbers 

< . . . < t_2 < t-i < to < ii < ^2 < • • • < 1 

with inf tfc = 0,suptfc = 1 such that 7o([i2J5i2i+i]) is a chord and 70 ([^2*- 1,^21]) belongs to Bd(C) 
for every i. 

If inf tfc = mintfc and supt^ = maxt^, then we have a finite sequence. In this case, we can 
approximate each arc 7o([t2j-i,i2j]) by a polygonal path 71 : [t2j-i,t2i] — ^ I^^ with the same 
endpoints, which avoids the interior of C. Letting 71 be identical to 70 on all shifting intervals 
[t2iit2i+i], we obtain a polygonal curve 71 consisting of finitely many line segments, and it can be 
easily achieved that 

^(7i)<^(7o) + e/2. 

Clearly, 71 is also a barrier for C, and it meets the requirements of the lemma. 

We can follow the same procedure when (tj) is a one-way or a two-way infinite sequence, except 
that now 71 consists of infinitely many segments. Suppose, for example, that there are infinitely 
many chords 7i([t2jii2j+i]) = lo{[t2i-,t2i+i])-, i = 0, 1,2, . . . . Obviously, these chords must converge 
to the common endpoint p of the curves 71,70, and 7. In this case, for a very large j, replace the 
final section 71 ([^2^, 1]) of 71 by a polygonal path connecting 7i(t2j) to 71(^2^+1), which avoids the 
interior of C. If j is large enough, then this transformation does not increase the length of the 
curve by more than e/A. Repeating the same trick at the other end of 71, if necessary, we obtain a 
polygonal (i. e., segment) barrier 72 for C, with 

%2)<^(7i) + 2(e/4)<%o)+e, 

as required. D 

Denote by per(C) the perimeter of a convex body C in the plane. The following lemma providing 
a lower bound on the length of an optimal barrier for C in terms of per(C), is used in the analysis 
of our approximation algorithms. Its proof is folklore; see e.g. |16j . 



Lemma 2. Let C he a convex body in the plane and let B he a harrier for C. Then the length of 
B is at least ^ ■ per(C). 

Proof. By Lemma [H we can assume w.l.o.g. that B is a segment barrier. Let B = {si, . . . ,Sn} 
consist of n segments of lengths ii = \si\, where L = \B\ = ^11=1^1- Let Oj G [0, vr) be the angle 
made by Sj with the x-axis. For each direction a G [0,7r), the blocking (opaqueness) condition for 
a convex body C can be written as 

n 

^^cos{a-ai)\>W{a). (1) 

i=l 

Here Ty(a) is the width of C in direction a, i. e., the minimum width of a strip of parallel lines 
enclosing C, whose lines are orthogonal to direction a. By integrating this inequality over the 
interval [0, vr], one gets: 



^ nil rTT 

y^^i \ cos{a — ai)\ da > / W{a) da. 
,_ 1 Jo Jo 



(2) 



According to Cauchy's surface area formula [36', pp. 283-284], for any planar convex body C, we 
have 

/ W{a) da = per(C). (3) 

Since 

/•TT 

/ I cos(a — ai)\ da = 2, 
Jo 

we get 

n 

2L = Y^ 2ii > per(C7) ^ L > - . per(C), (4) 

as required. D 

For instance, for the square, per(C) = 4, and Lemma [2] immediately gives L > 2, the lower 
bound of Jones |24j ) . 

Remark. Obviously, the boundary of C, dC, is a barrier for C of length per(C). Consequently, 
once Lemma [2] is established, a 2-approximation (for each type of barrier) follows immediately. 

A key fact in the analysis of our approximation algorithms is the following lemma. This in- 
equality is implicit in |44] : another proof can be found in [12J. 

Lemma 3. Let P be a convex polygon. Then the minimum-perimeter rectangle R containing P 
satisfies per(i?) < ;^per(P). 

Let P be a convex polygon with n vertices. Let OPTarb(-P)) OPTconn(P) and OPTarc(-P) denote 
optimal barrier lengths of the types arbitrary, connected, and single-arc. Observe the following 
inequalities: 

OPTa,b(i^) < OPTconn(i^) < OPTarc(i^). (5) 

We first deal with connected barriers, and then with arbitrary (i. e., possibly disconnected) 
barriers. 

3 Connected barriers 

Theorem 1. Given a convex polygon P with n vertices, a connected polygonal harrier whose length 
is at most 1.5716 times longer than the optimal can he computed in 0{n) time. 

Proof. We start with the following algorithm Al that computes a connected barrier consisting 
of a single-arc; refer to Fig. [21 First compute a parallel strip of minimum width enclosing P. 
Assume w.l.o.g. that the strip is bounded by the two horizontal lines ii and £2- Second, compute 
a minimal orthogonal {i. e., vertical) strip enclosing P, bounded by the two vertical lines ^3 and 
£4. Let a, 6, c, d, e, / be the six segments on £3 and £4 as shown in the figure; here b and e are the 
two (possibly degenerate) segments on the boundary of P. Let Pi be the polygonal path (on P's 
boundary) between the lower vertices of b and e. Let P2 be the polygonal path (on P's boundary) 
between the top vertices of b and e. 

Consider the following two barriers for P: Pi consists of the polygonal path Pi extended upward 
at both ends until they reach £2. P2 consists of the polygonal path P2 extended downwards at both 
ends until they reach ii. The algorithm returns the shorter of the two. We show below that its 
approximation ratio is at most ^^^ = 1.5834. . .. 
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Figure 2: The approximation algorithm Al returns Bi (in bold lines). 



Let p, u;, and r, respectively, be the perimeter, the width, and the in-radius of P . Clearly 



|i^i| + |-P2| + |6| + |e| 



V- 



We have the following equalities: 



|Bi| = |a| + |6| + |Pi| + |e| + |/|, 
I-B2I = |c| + |6| + |P2| + |e| + |d|. 

By adding them up we get 

l^il + IB2I = |Pi| + IP2I + \h\ + |e| + 2u; = p + 2w. 

Hence 

min{|Bi|,|52|} <p/2 + u;. 

By Blaschke's Theorem [6] (see also [451 Exercise 2-5]), every planar convex body of width w 
contains a disk of radius w/3, hence r > wj'i. This inequality cannot be improved: equality is 
attained for the equilateral triangle. According to a result of Eggleston [13], the optimal connected 
barrier for a disk of radius r has length (7r + 2)r. It follows that the optimal connected barrier for P 
has length at least (vr + 2)t(;/3. By Lemma [21 p/2 is another lower bound on the optimal solution. 
Thus the approximation ratio of the algorithm Al is at most 



p/2 + w 



max{(7r + 2)t(;/3, p/2} 



mm 



mm 



p/2 + w p/2 + w 



(7r + 2)w/3' p/2 
3 



2(7r + 2) 



p 3 

w vr + 2 



1 + 2 



w 



P 



One can check that the quadratic equation 



3x 



+ 



2{tt + 2) 7r + 2 



1 + 



has one positive real root 



xo 



2(7r + 2) 



Consequently, the approximation ratio of the algorithm Al is at most 1 + ■:^^-r2 — f+2 ~ 
1.5834.... Clearly the algorithm takes 0{n) time, since computing the width of P takes 0{n) 
time [371142], and the two barriers Bi and B2 can be computed within the same time. 
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We next achieve a better approximation, 1.5716, by means of a more elaborated approach. The 
idea is to do something different when P is "close to" an equilateral triangle. In this case, one of 
the two barriers Bi and B2 computed by algorithm Al is substantially shorter than the average 
of the two, namely, min{|Si|, I-B2I} is substantially shorter than (|-Bi| + |i?2|)/2, and the previous 
argument becomes wasteful. Our revised algorithm is A2. 

To explain the algorithm, we need to enter the details of the proof of Blaschke's Theorem, as 
given in |45^ Exercise 6-2]. Let Q be a largest circle contained in P; let r be its radius. Then Q 
either contains two diametrically opposite points, or else it contains three boundary points of P 
which form an acute-angled triangle. In the former (easier) case, r = w/2, and this yields a much 
better approximation than that obtained earlier using the inequality r > w/3; to put it short, this is 
not the bottleneck case. Assume therefore that Q is incident to three boundary points, A,B,CgP 
which form an acute triangle, AABC. Then the supporting lines at A, B, C must form a triangle 
T = AA'B'C which is circumscribed to both the polygon P and the circle O. Denote the sides 
of this triangle by o, 6, c, where o is a longest side, and the corresponding altitudes by ha, hf,, he- 
Denote by Wa, Wb, Wc the widths of P in the directions of a, 6 and c, respectively. Obviously, we 
have ha > Wa > w, hi, > w^ > w, and hc>Wc> w. See Fig. [3l 




Figure 3: P (in bold lines) and T. 



We now present the revised algorithm. The algorithm A2 first computes the two barriers Bi 
and B2 as done by algorithm Al. In addition, it also computes a third barrier, B^, which is a 
Steiner minimal tree of the three points A' , B' , C", if they exist; otherwise -B3 is undefined and 
|i?3| = 00. Since P is contained in T, B3, which is a connected barrier for T, is also a connected 
barrier for P. The algorithm then returns the shorter of the three barriers, Bi, B2, -B3. 

Recall that a Steiner minimal tree of three points that determine no angle larger or equal to 
27r/3 is a star, whose any two consecutive edges make an angle of 27r/3 between them; see e. g., |21j . 
or [Ml Ch. 6]. 

Returning now to the proof of Blaschke's Theorem, the area of T can written is several ways: 

{a + b + c)r aha bh^ chc 
2 ~ ~Y' ~~ ~ ~' 

Since a > 5, a > c it follows that 



Area(r) 



(6) 



h„ 



w„ 



.ha>^>^>^ 



(7) 
a + b + c "^ 3 3 3 

This concludes the proof of Blaschke's Theorem. Observe that if a is somewhat larger than 
{a + b + c)/3, then r is somewhat larger than w/3, and one could use this improved bound to get a 
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better approximation ratio as in the analysis of Case 1. We next analyze the approximation ratio 
of algorithm A2. We can assume w.l.o.g. that the perimeter of T is 1, i. e., a + b + c = 1, and 
further that c < b. Then a > 1/3. We make use of two parameters A and 6 = 3 — 1/A, where 
1/3 < A < 3/8 and correspondingly < 5 < 1/3, which will be later set to A = 0.3403 . . . and 
6 = 0.0615 ... in order to optimize the approximation ratio of A2. We distinguish two cases: 

Case 1. a > X. Then according to ([7]), we have r = aha ^ o,Wa > aw > Xw. Similar to the 
previous analysis of Al, the approximation ratio of A2 is at most 

p/2 + w ( p/2 + w p/2 + w 

— min 



max{(7r + 2)At(;, p/2} [{tt + 2)Xw' p/2 



1 p 1 w 

mm < — — 7 • h -— -, 1 + 2 • — 

\2A(7r + 2) w A(7r + 2)' p 



As before, one can easily check that the quadratic equation 

X 1 2 

+ . . . ^^ = 1 + - 



2A(7r + 2) A(7r + 2) x 

has one positive real root 

xo = 2A(7r + 2). 

Under the assumption in Case 1, it follows that the approximation ratio is at most 

2 1 

1 + — = 1 + 



xo A(7r + 2) 

For future reference, set 

Case 2. a < X. Obviously, we also have b,c < X. Then b = l — a — c>l — 2A, and similarly, 
c > 1 — 2A. To summarize, 

1-2A < a,6,c< A, a > -. (9) 

Recall that a<A<3/8<\/2-l, which implies that a^ < 2(1 - a)2/4 < 6^ + c^. It follows that 
T is an acute triangle. We further distinguish two sub-cases. Case 2.1 and Case 2.2. 

Case 2.1. At least one of the following three inequalities holds: (i) Wa < {I — 6)ha; (ii) 
Wb < {I — S)hb; (iii) Wc < {I — 6)hc. Let ^ £ {a,b,c} and assume that w^ < {1 — 5)h^. Then ([6]) 
and (jlD yield 

r = (h^>{l- 2A) Y^ > (1 - 2A) jj^ = Xw. 

As in the analysis of Case 1, it follows that the approximation ratio is again at most pi under the 
assumption in Case 2.1. 

Case 2.2. None of the inequalities in Case 2.1 holds. We then have w^ > {1 — S)h^, for each 
^ € {a,b,c}. Construct a triangle T' containing P as described below, and as shown in Fig. SI 
Assume w.l.o.g. that the side a is the horizontal base of T. Consider the three lines, ia, ^b, ^c, each 
parallel to the corresponding side of T: a line ia parallel to a and tangent to P from above, etc. 
Observe that the triangles T and T' are similar, by construction. Let 6ia, 52b, 63c, be the segments 
of intersection of the three lines with T. 



-r 



Figure 4: P (in bold lines), T (in solid lines) and T' (in dashed lines). 

By the assumption of Case 2.2., we have h'c < 5h^, for each ^ G {a, 6, c}. where h'^, h'j^, h'^ denote 

the altitudes from A' , B' , and C" in the three smaller similar triangles incident to A', B' , and C". 

It follows that 

h' h' h' 

ll'r, I'-h ll'^ 

5l + 02 + 03 



5i + d2 + 63='4^ + 'T^+'-^<36. 



ha hh he 

It is easily seen that the similarity ratio between T' and T is 2 — 5i — 62 — S3. By the previous 
bound, this ratio is at least 2 — 36 > 1. Observe that P is incident to the three sides of the acute 
triangle T' . It is well-known that the minimum-perimeter triangle inscribed in a given acute triangle 
A (i. e., with a vertex incident to each side of A) is the orthic triangle of A |271 Theorem 17]; or 
see [391 Ch. 5]. The vertices of the orthic triangle are the feet of the altitudes of A. It is also known 
that the semiperimeter of the orthic triangle of an acute triangle with semiperimeter s, and sides 
X, y and z is equal to 

4s(s — x){s — y){s — z) 

xyz 
In particular, since a + h + c= 1, the semiperimeter of the orthic triangle of T is 

2{\-a){\-h){\-c) ^ 
abc 

Since the similarity ratio between T' and T is at least 2 — 35, by taking into account Q, we obtain 
that the semiperimeter of the orthic triangle of T' is at least 

2(2-35)(l-a)(i-6)(i-c) ^(3_A / J_ _ \ / J_ _ \ /J_ _ ^ 



abc VA / \2a J \2b J \2c 

3 J\ / I ^ ^ 



Since P is incident to the three sides of T' , its semiperimeter p/2 is bounded from below by the 
above expression, thus 
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We now bound from above the length of the third barrier B^. Recall that a > b > c. We have 
^C < 7r/3, thus -cos(ZC + 7r/3) < 1/2. One can deduce from [in Section 5] (or from [391 Ch. 6]) 
and by using our assumptions in Case 2 that 

|S3|2 = a2 + 62_2a6cos(ZC + 7r/3) <a'^ + b'^ + ab< 3A^ 

hence |i?3| < Av3- Taking into account (jlOp . under the assumptions in Case 2.2, the approximation 
ratio is at most 

^2 := —yz TTl Ts- (11) 

Clearly, the approximation ratio of algorithm A2 is at most p = maxjpi, /32}. To balance Cases 
1 and 2.1 with Case 2.2, we let A be the solution to the equation pi{X) = /02(A) below; recall ([8]) 
and ([n]): 

l + -i-^ = '-^ ,. (12) 

A routine calculation shows that A = 0.3403 . . . and, correspondingly, 5 = 3 — 1/A = 0.0615 . . . and 
Pi = P2 = 1-5715 .... We conclude that the approximation ratio of algorithm A2 is at most 1.5716, 
as claimed. 

The largest circle inscribed in a convex polygon can be found by linear programming in linear 
time [35]. Computing By, given T takes constant time, thus B^ can be computed in 0{n) time. 
Recall that Bi and B2 can be computed in 0{n) time too. Consequently, the algorithm A2 takes 
0{n) time. D 

It is easy to see that the connected barrier computed by A2 is not optimal in general (in the 
class of connected barriers). The square gives an easy example. The length of the third barrier from 
the left in Fig. [T]is 1 + -v/3, while the length of the barrier computed by A2 is 3 {\Bi\ = \B2\ = 3, 
\B^\ = 00). This example shows a lower bound of 1.098... on the approximation ratio of the 
algorithm A2. 

4 Single-arc barriers 



Since algorithm Al computes a single-arc barrier, and we have OPTconn(^) < OPTarc(-P)) we 

~2 



immediately get an approximation algorithm with ratio ^^^^ = 1.5834. . . for computing single- arc 



barriers. 

Theorem 2. Given a convex polygon P with n vertices, a single-arc polygonal barrier whose length 
is at most ^i| = 1.5834. . . times longer than the optimal can be computed in 0{n) time. 

One may ask whether the single arc barrier computed by Al is optimal (in the class of single arc 
barriers). We show that this is not the case: Consider (a sufficiently fine polygonal approximation 
of) a Reuleaux triangle T of (constant) width 1, with three vertices a, 6, c. Now slightly shave 
the two corners at b and c and obtain a convex body T' of (minimum) width 1 — e along be. The 
algorithm Al would return a curve of length close to 7r/2 + 1 = 2.57. . ., while the optimal curve 
has length at most 27r/3 + 2 — ^/3 = 2.36 .... This example shows a lower bound of 1.088 ... on the 
approximation ratio of the algorithm Al. On the other hand, we believe that the approximation 
ratio of Al is much closer to this lower bound than to 1.5834 .... 

We next present an improved version A3 of our algorithm Al that computes the shortest 
single-arc barrier of the form shown in Fig. [2l Let P be a convex polygon with n sides, and let i 
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be a line tangent to the polygon, i. e., P Di consists of a vertex of P or a side of P. For simplicity 
assume that i is the x-axis, and that P lies in the closed halfplane y >0 above i. Let T = (^1,^2) 
be a minimal vertical strip enclosing P. Let pi € iiCiP and p2 G ^2 H P, be the two points of P of 
minimum y-coordinates on the two vertical lines defining the strip. Let qi £ ii and (72 £ ^2 be the 
projections of pi and p2, respectively, on i, and avc{pi,p2) C dP be the polygonal arc connecting 
pi and P2 on the top boundary of P . 

The [/-curve corresponding to P and £, denoted U{P,i) is the polygonal curve obtained by 
concatenating qipi, arc{pi,p2), andp292) in this order. Obviously, for any line £, the curve U{P,i) 
is a single-arc barrier for P. Let [/min(-P) be the [/-curve of minimum length over all directions 
a G [0,7r) {i. e., lines £ of direction a). 

We next show that given P, the curve [/min(-P) can be computed in 0{n) time. The algorithm 
A3 is very simple: instead of rotating a line I around P, we fix £ to be horizontal, and rotate P over 
£ by one full rotation (of angle 27r). We only compute the lengths of the [/-curves corresponding to 
lines i, £1, £2, supporting one edge of the polygon. The [/-curve of minimum length among these 
is output. There are at most 3n such discrete angles (directions), and the length of a [/-curve for 
one such angle can be computed in constant time from the length of the [/-curve for the previous 
angle. The algorithm is similar to the classic rotating calipers algorithm of Toussaint [42], and it 
takes 0{n) time by the previous observation. 

To justify its correctness, it suffices to show that if each of the lines £, £1, £2 is incident to only 
one vertex of P, then the corresponding [/-curve is not minimal. 

Lemma 4. Let P be a convex polygon tangent to a line £ at a vertex v £ P only, and tangent to 
£1 and £2 at vertices pi and p2 only. Then the corresponding U -curve U{P, £) is not minimal. 

Proof. For convenience, assume that £ is horizontal, and that P lies in the closed halfplane above 
£. Refer to Fig. El 




qi V '72 

Figure 5: The curve U{P,£). 

Let pi,qi G £\ and ^2,^2 £ £2 be as defined earlier. Observe that v belongs to the closed 
segment q\q2. If w = gi or f = 52, then by rotating P (clockwise or counterclockwise, as needed) 
around i; by a small angle, the length of the curve U{P,£) decreases. So we can assume that v 
lies in the interior of the segment qiq2. Observe that if P rotates clockwise or counterclockwise by 
a small angle around v, pi and p2 remain the same, so the angle Zpivp2 also remains the same. 
Put a = Zqivpi, j3 = Zpii;p2, and 7 = Zp2vq2, so a + /? -|- 7 = vr. Put a = \piv\, b = \piP2\, and 
c = \vp2\. The length of U{P,£) for this angle a is 

/(a) = asina -|- |arc(pi,p2)| + csin7. 
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The first two derivatives of /(•) are 

f'{a) = a cos a — ccos{-k — a — /3) = a cos a — CC0S7. 
f '{a) = —a sin a — csin(7r — a — (3) = —a sin a — csin7. 

Since 0,7 S (0, 7r/2), we have /"(a) < 0, wliich means that f{a) is not a local minimum. D 

Theorem 3. Given a convex polygon P with n vertices, the single- arc barrier (polygonal curve) 
Umin{P) can be computed in 0{n) time. 

Obviously, the approximation ratio of the algorithm A3 is not worse than that achieved by 
algorithm Al, hence it is also bounded by ^x| = 1-5834 .... One may ask again whether the single 
arc barrier computed by A3 is optimal (in the class of single arc barriers). We show again that this 
is not the case. Consider the pentagon with vertices (0,e), (—3,0), (— 1,— e), (1,— e), (3,0). The 
optimal curve is no longer than the curve ((—3,0), (—1, — e), (0,e), (1, — e), (3,0)), whose length is 
6 + O(e^). On the other hand, the algorithm A3 returns a curve of length 6 + r2(e). Let now e be 
sufficiently small. 

We can fine-tune (numerically) the above pentagon to obtain a lower bound of 1.065 ... on the 
approximation ratio of A3. See Fig. [6j 




(a) 




(b) 




(c) 



Figure 6: A pentagon with five vertices A = (0,/i), B = {—x,y), C = (—1,0), D = (1,0), E = {x,y), where 
X = 1.4507 ..., y = 0.2072 . . ., and /i = 0.3806 . . .. The algorithm A3 returns a barrier of length 3.3364 . . . 
as in (a) or (b), but the barrier in (c) has a shorter length of 3.132 . . .. This gives a lower bound of 1.065 . . . 
on the approximation ratio of the algorithm. 



5 Arbitrary barriers 

Theorem 4. Given a convex polygon P with n vertices, a (possibly disconnected) segment barrier 
for P , whose length is at most 2 + ~^^ = 1.5867 . . . times longer than the optimal can be computed 
in 0{n) time. 

Proof. Consider the following algorithm A4 which computes a (generally disconnected) barrier. 
First compute a minimum-perimeter rectangle R containing P; refer to Fig. [71 Let a,b,c,d,e,f,g,h, 
i,j,k,l be the 12 segments on the boundary of R as shown in the figure; here 6, e, h and k are 
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(possibly degenerate) segments on the boundary of P contained in the left, bottom, right and top 
side of R. Let Pi, z = 1, 2, 3, 4 be the four polygonal paths on P's boundary, connecting these four 
segments as shown in the figure. 

D 





A d 



e f B A 

Figure 7: The approximation algorithm A4. 



B 



Consider four barriers for P, denoted Bi, for i = 1,2,3,4. Bi consists of the polygonal path 
Pi extended at both ends on the corresponding rectangle sides, and the height from the opposite 
rectangle vertex in the complementary right angled triangle; see Fig. [Tl^right). The algorithm 
returns the shortest of the four barriers. Let h^, hs, he, ho denote the four altitudes from A, B, 
C, and D, respectively, in the right-angled triangles /lABD, ABCA, ACBD, and ABAC. We 
have |/ia| = |/ib| = \hc\ = {hoi and the following other equalities: 

|Bi| = |a| + |6| + |Pi| + |e| + |/| + |/ic|, 
\B2\ = \d\ + \e\ + IP2I + \h\ + |i| + {hnl, 

\B3\ = \g\ + \h\ + \P3\ + \k\ + \l\ + \hA\, 

\B^\ = \j\ + \k\ + \P4\ + \b\ + \c\ + \hB\. 

By adding them up yields 

4 4 

Y, \B^\ = {\b\ + |e| + \h\ + \k\ + Y, \P^\) + (|a| + • • • + |/|) + (|/ia| + Ihsl + \hc\ + Ihnl) 

= per(P)+per(i?) + 4|/iA|. (13) 

The length of the altitude |/ia| in the right-angled triangle AABD is given by the formula 

\hA\ 



xy 



where x and y are the lengths of the two sides of R. By Lemma [3] we have 

4 
per(it') = 2(x + 2/) < - per(P). 



Under this constraint, |/ia| is maximized for x 

per(P) 



y 



vr 

per(P) 



namely 



\hA\< 



^T^/2 



4|/iA|<^per(P). 



■K 



Hence ()T3|) yields 



1 / A /9 \ 

minlSJ < - IH h -^ per(P). 

i 4 V vr vr / 
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Recall that per(P)/2 is a lower bound on the weight of an optimal solution. The ratio between 
the length of the solution and the lower bound on the optimal solution is 

. + 4 + 2^2^1 2W2^ 

27r 2 vr 

Consequently, the approximation ratio of the algorithm A4 is 2 + ^^ = 1.5867 .... The algorithm 
takes 0{n) time, since computing the minimum-perimeter rectangle containing P takes 0{n) time 
with the standard technique of rotating calipers |371 [52] . This completes the proof of Theorem^ D 

The above analysis of the approximation ratio of A4 is tight for a circle, in particular if P is a 
regular n-gon with n tending to infinity. Indeed, for a unit-radius circle, per(P)/2 = tt while the 
length of the barrier computed by A4 is 2 + 7r/2 + v2. 

6 Interior-restricted versus unrestricted barriers 

In certain instances, it is infeasible to construct barriers guarding a specific domain outside the do- 
main (which presumably belongs to someone else). We call such barriers constrained to the interior 
and the boundary of the domain, interior-restricted, or just interior, and all others unrestricted. 
For example, all four barriers for the unit square illustrated in Fig. [T] are interior barriers. 

In the late 1980s, Akman [T] soon followed by Dublish [11] had reported algorithms for com- 
puting a minimum interior-restricted barrier of a given convex polygon (they refer to such a barrier 
as an opaque minimal forest of the polygon). Both algorithms however have been shown to be 
incorrect by Shermer pij in 1991. He also proposed (conjectured) a new exact algorithm instead, 
but apparently, so far no one succeeded to prove its correctness. To the best of our knowledge, the 
computational complexity of computing a shortest barrier (either interior-restricted or unrestricted) 
for a given convex polygon remains open. 

Next we show that a minimum connected interior barrier for a convex polygon can be computed 
efficiently: 

Theorem 5. Given a convex polygon P, a minimum Steiner tree of the vertices of P forms a mini- 
mum connected interior harrier for P . Consequently, there is a fully polynomial-time approximation 
scheme for finding a minimum connected interior harrier for a convex polygon. 

Proof. Let B be an optimal barrier. For each vertex v G P, consider a line i^ tangent to P at v, 
such that P Ciiy = {v}. Since B lies in P, (.^ can be only blocked by w, so f G B. Now since B is 
connected and includes all vertices of P, its length is at least that of a minimum Steiner tree of P, as 
claimed. Recall that the minimum Steiner tree problem for n points in the plane in convex position 
admits a fully polynomial-time approximation scheme that achieves an approximation ratio of 1 + e 
and runs in time 0{n^ /e^) for any e > [38]. D 

A minimum single-arc interior barrier for a convex polygon can be also computed efficiently. 
As it turns out, this problem is equivalent to that of finding a shortest traveling salesman path 
(i. e., Hamiltonian path) for the n vertices of the polygon. 

Theorem 6. Given a convex polygon P, a minimum Hamiltonian path of the vertices of P forms 
a minimum single-arc interior harrier for P. Consequently, there is an 0(n'^)-time exact algorithm 
for finding a minimum single-arc interior harrier for a convex polygon with n vertices. 
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Proof. The same argument as in the proof of Theorem [5] shows that any interior barrier for P 
must include all vertices of P. By the triangle inequality, the optimal single-arc barrier visits each 
vertex exactly once. Thus a minimum Hamiltonian path of the vertices forms a minimum single-arc 
interior barrier. 

We now present a dynamic programming algorithm for finding a minimum Hamiltonian path 
of the vertices of a convex polygon. Let {vq, . . . ,Vn-i} be the n vertices of the convex polygon 
in counter-clockwise order; for convenience, the indices are modulo n, e.g., Vn = vq. Denote by 
dist(i, j) the Euclidean distance between the two vertices Vi and Vj. For the subset of vertices from 
Vi to Vj counter-clockwise along the polygon, denote by S{i, j) the minimum length of a Hamiltonian 
path starting at Vi, and denote by T(i,j) the minimum length of a Hamiltonian path starting at 
Vj. Note that a minimum Hamiltonian path must not intersect itself. Thus the two tables S and 
T can be computed by dynamic programming with the base cases 

S{i, i + 1) = T{i, i + l)= dist(«, i + 1) 

and with the recurrences 

S{i,j) = min{dist(z,i + 1) + S{i + l,j), dist{i,j) +T{i + l,j)}, 
T(i,j) = min{dist(j,j - 1) +T{i,j - 1), dist(j,z) + S{i,j - 1)}. 

Then the minimum length of a Hamiltonian path on the n vertices is 

minmin{dist(«,i + 1) + S{i + l,z — 1), dist{i,i — 1) + T{i + l,i — 1)}. 

i 

The running time of the algorithm is clearly 0{v?). D 

Remark. Observe that the unit square contains a disk of radius 1/2. According to the result of 
Eggleston mentioned earlier [13], the optimal (not necessarily interior-restricted) connected barrier 
for a disk of radius r has length (vr -|- 2)r. This optimal barrier is a single curve consisting of 
half the disk perimeter and two segments of length equal to the disk radius. It follows that the 
optimal (not necessarily interior-restricted) connected barrier for the unit square has length at least 
(vr + 2)/2 = 7r/2 -|- 1 = 2.5707. . .. Compare this with the current best construction (illustrated in 
Fig. [H third from the left) of length 1 -|- yS = 2.7320 .... Note that this third construction in Fig. [1] 
gives the optimal connected interior barrier for the square because of Theorem [5l Further note that 
the first construction in Fig. [1] gives the optimal single-arc interior barrier because of Theorem [6l 

7 Conclusion 

Interesting questions remain open regarding the structure of optimal barriers and the computational 
complexity of computing such barriers. For instance: 

(1) Does there exist an absolute constant c > (perhaps zero) such that the following holds? 
The shortest barrier for any convex polygon with n vertices is a segment barrier consisting of 
at most n + c segments. 

(2) Is there a polynomial-time algorithm for computing a shortest barrier for a given convex 
polygon with n vertices? 

(3) Can one give a characterization of the class of convex polygons whose optimal barriers are 
interior? 
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In connection with question (2) above, let us notice that the problem of deciding whether a 
given segment barrier B is an opaque set for a given convex polygon is solvable in polynomial time: 

Theorem 7. Given a convex polygon P with n vertices, and a segment barrier B with k segments, 
there is a polynomial-time algorithm for deciding whether B is an opaque set for P. 

Proof. Let V{B) denote the 2k endpoints of the segments in B. Consider the set of lines (directions) 
C determined either by pairs of distinct points in V{B) or that are incident to a point in V{B) and 
tangent to P. Observe that C has 0{k'^) elements, and it can be easily constructed in 0{nk + fc^) 
time. If B is not an opaque set for P, there exists a line i £ C such that a small rotation (clockwise 
or counterclockwise) around i yields a direction, say i'^ or i~ , such that the projection of B onto 
the line orthogonal to it does not cover the projection of P onto the same line. That is, the union 
of the projection segments does not include the segment which represents the projection of P. We 
say that the opaqueness condition fails with respect to £"*" or i~ . 

To see this, take a line that intersects P without intersecting B. Fix a point p in P (in the 
interior or on the boundary of P) on this line and rotate the line around p until it hits a segment 
in B, say at its endpoint q. Start rotating the line around q until either it becomes tangent to P 
(as it leaves P), or it hits another segment endpoint in V{B). 

For a given i £ C the opaqueness condition for £+ and i~ can be easily checked in 0{n + k) time. 
Since there are 0{k'^) lines in C, the overall opaqueness can be checked in 0{{n+k)-k'^) time. Hence 
whether B is an opaque set for P can be determined in 0{nk + k'^ + {n + k) ■ k"^) = 0{{n + k) ■ k'^) 
time. (A faster algorithm can be obtained by using rotational sweep [H p. 328].) D 

We have presented several approximation and exact algorithms for computing shortest barriers 
of various kinds, for a given convex polygon. The two approximation algorithms with ratios close 
to 1.58 probably cannot be improved substantially without either increasing their computational 
complexity or finding a better lower bound on the optimal solution than that given by Lemma [21 
The question of finding a better lower bound is particularly intriguing, since even for the simplest 
polygons, such as a square, we don't possess any better tool. While much research up to date focused 
on upper or lower bounds for specific example shapes, obtaining a polynomial time approximation 
scheme (in the class of arbitrary barriers) for an arbitrary convex polygon is perhaps not out of 
reach. 

References 

[1] V. Akman, An algorithm for determining an opaque minimal forest of a convex polygon, 
Information Processing Letters, 24 (1987), 193-198. 

[2] D. Asimov and J. L. Gerver, Minimum opaque manifolds, Geom. Dedicata 133 (2008), 67-82. 

[3] F. Bagemihl, Some opaque subsets of a square, Michigan Math. J. 6 (1959), 99-103. 

[4] I. Barany and Z. Fiiredi, Covering all secants of a square, in Intuitive Geometry (G. Fejes 
Toth, editor), Colloq. Math. Soc. Janos Bolyai, vol. 48 (Siofok, Hungary, 1985), pp. 19-27, 
North-Holland, Amsterdam, 1987. 

[5] M. de Berg, O. Cheong, M. van Kreveld, and M. Overmars, Computational Geometry, 3rd 
edition. Springer Verlag, 2008. 

[6] W. Blaschke, Kreis und Kugel, 2te Aufl., Walter de Gruyter, Berlin, 1956. 

16 



[7] K. A. Brakke, The opaque cube problem, American Mathematical Monthly, 99(9) (1992), 
866-871. 

[8] H. T. Croft, Curves intersecting certain sets of great-circles on the sphere, J. London Math. 
Soc. (2) 1 (1969), 461-469. 

[9] H. T. Croft, K. J. Falconer, and R. K. Guy, Unsolved Problems in Geometry, Springer, New 
York, 1991. 

[10] E. D. Demaine and J. O'Rourke, Open problems from CCCG 2007, in Proceedings of the 20th 
Canadian Conference on Computational Geometry (CCCG 2008), Montreal, Canada, August 
2008, pp. 183-190. 

[11] P. Dublish, An 0{n^) algorithm for finding the minimal opaque forest of a convex polygon, 
Information Processing letters, 29(5) (1988), 275-276. 

[12] A. Dumitrescu and M. Jiang, Minimum-perimeter intersecting polygons, Algorithmica, to ap- 
pear. An extended abstract in Proceedings of the 9th Latin American Theoretical Informatics 
Symposium (LATIN 2010), Oaxaca, Mexico, April 2010, LNCS Vol. 6034, pp. 433-445. 

[13] H. G. Eggleston, The maximal in-radius of the convex cover of a plane connected set of given 
length, Proc. London Math. Soc. (3) 45 (1982), 456-478. 

[14] P. Erdos and J. Pach, On a problem of L. Fejes Toth, Discrete Mathematics, 30(2) (1980), 
103-109. 

[15] V. Faber and J. Mycielski, The shortest curve that meets all the lines that meet a convex 
body, American Mathematical Monthly, 93 (1986), 796-801. 

[16] V. Faber, J. Mycielski and P. Pedersen, On the shortest curve that meets all the lines which 
meet a circle, Ann. Polon. Math., 44 (1984), 249-266. 

[17] L. Fejes Toth, Exploring a planet, American Mathematical Monthly, 80 (1973), 1043-1044. 

[18] L. Fejes Toth, Remarks on a dual of Tarski's plank problem. Mat. Lapok., 25 (1974), 13-20. 

[19] S. R. Finch, Mathematical Constants, Cambridge University Press, 2003. 

[20] M. Gardner, The opaque cube problem, Cubism for Fun 23 (March 1990), p. 15. 

[21] E. N. Gilbert and H. O. Pollak, Steiner minimal trees, SIAM Journal on Applied Mathematics, 
16 (1968), 1-29. 

[22] H. M. S. Gupta and N. C. B. Mazumdar, A note on certain plane sets of points, Bull. Calcutta 
Math. Soc. 47 (1955), 199-201. 

[23] R. Honsberger, Mathematical Morsels, Dolciani Mathematical Expositions, No. 3, The Math- 
ematical Association of America, 1978. 

[24] R. E. D. Jones, Opaque sets of degree a, American Mathematical Monthly, 71 (1964), 535-537. 

[25] H. Joris, Le chasseur perdu dans le foret: une probleme de geometrie plane, Elemente der 
Mathematik, 35 (1980), 1-14. 



17 



[26] B. Kawohl, Some nonconvex shape optimization problems, in Optimal Shape Design (A. Cellina 
and A. Ornelas, editors), vol. 1740/2000 of Lecture Notes in Mathematics, Springer, 2000. 

[27] N. D. Kazarinoff, Geometric Inequalities, Random House, New York, 1961. 

[28] W. Kern and A. Wanka, On a problem about covering lines by squares. Discrete and Compu- 
tational Geometry, 5 (1990), 77-82. 

[29] R. Klotzler, Universale Rettungskurven I, Zeitschrifte fiir Analysis und ihre Anwendungen, 5 
(1986), 27-38. 

[30] R. Klotzler and S. Pickenhain, Universale Rettungskurven II, Zeitschrifte fiir Analysis und 
ihre Anwendungen, 6 (1987), 363-369. 

[31] E. Kranakis, D. Krizanc, L. Narayanan, K. Xu, Inapproximability of the perimeter defense 
problem, in Proceedings of the 21st Canadian Conference on Computational Geometry (CCCG 
2009), Vancouver, Canada, August 2009, pp. 153-156. 

[32] E. Makai, Jr., On a dual of Tarski's plank problem. Discrete Geometric, 2, Kolloq., Inst. Math. 
Univ. Salzburg, 1980, pp. 127-132. 

[33] E. Makai, Jr. and J. Pach, Controlling function classes and covering Euclidean space, Studia 
Scientiarum Mathematicum Hungaricae, 18 (1983), 435-459. 

[34] S. Mazurkiewicz, Sur un ensemble ferme, punctiforme, qui rencontre toute droite passant par 
un certain domaine (Polish, French summary), Prace Mat.-Fiz. 27 (1916), 11-16. 

[35] N. Megiddo, Linear programming in linear time when the dimension is fixed. Journal of the 
ACM, 31 (1984), 114-127. 

[36] J. Pach and P. K. Agarwal, Combinatorial Geometry, John Wiley, New York, 1995. 

[37] F. Preparata and M. I. Shamos, Computational Geometry, Springer- Verlag, New York, 1985. 

[38] J. S. Provan, Convexity and the Steiner tree problem. Networks, 18 (1988), 55-72. 

[39] H. Rademacher and O. Toeplitz, The Enjoyment of Mathematics, Princeton University Press, 
1957. 

[40] T. Richardson and L. Shepp, The "point" goalie problem. Discrete and Computational Geom- 
etry, 20 (2003), 649-669. 

[41] T. Shermer, A counterexample to the algorithms for determining opaque minimal forests. 
Information Processing letters, 40 (1991), 41-42. 

[42] G. Toussaint, Solving geometric problems with the rotating calipers. Proceedings of Mediter- 
ranean Electrotechnical Conference (MELECON '83), Athens. 

[43] P. Valtr, Unit squares intersecting all secants of a square, Discrete and Computational Geom- 
etry, 11 (1994), 235-239. 

[44] E. Welzl, The smallest rectangle enclosing a closed curve of length vr, manuscript, 1993. Avail- 
able at http : //www . inf . ethz . ch/personal/emo/SmallPieces . html. 

[45] I. M. Yaglom and V. G. Boltyanskii, Convex Figures, Holt, Rinehart and Winston, New York, 
1961. 



18 



